Lateral electron localization by the induced surface charge 
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We investigate the problem of the electron interacting with the charge induced on the metal or 
dielectric surface. We show that the interaction between the electron and the induced surface charge 
leads to the lateral confinement of the electron. As a result the electron propagates parallel to the 
surface not as a plane wave but as a wave packet of a Gaussian shape. The electron moving together 
with the induced charge can be treated as a new quasi-particle, which we call inducton. We discuss 
a possible experimental evidence for inductons in semiconductor nanostructures, metal- vacuum, and 
dielectric-vacuum interfaces. 
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The electron interacting with the charge induced on 
the metal (dielectric) surface was a subject of many ex- 
perimental H S H H H and theoretical 0, H ll ES 
llll papers. In semiconductor nanostructures, the 
induced potential leads to a screening of the electron- 
electron interaction and a modification of the potential 
confining the charge carriers in the quantum wells, wires, 
and dots HI El El El, E3 • In papers 0,11 the classical 
image-charge potential was applied to study the electron 
localization in the direction normal to the metal surface. 
The lateral localization of the electron, i.e., parallel to 
the metal surface, has been observed by the time-resolved 
two-photon photoemission [j, ^bl . 

In the present Letter we study the electron moving 
near the surface of the metal or dielectric. We show that 
the interaction with the induced charge leads to the lat- 
eral confinement potential acting on the electron. In the 
existing literature, e.g., Ill la, Efll 111 ll3 , this potential 
has been overlooked due to the neglect of the dependence 
of the induced charge distribution on the electron wave 
function. 

We first formulate a two-dimensional (2D) model for a 
typical metal-semiconductor planar structure. Next, we 
extend the model to the three-dimensional (3D) system 
and compare the results of calculations with the available 
experimental data. We show that the observed 0, 0] lat- 
eral electron localization can be explained by the effect 
of the induced surface charge without a necessity of tak- 
ing into account neither the surface corrugation nor the 
electron-dipol coupling. 

We consider a planar structure, in which an electron 
is confined in the quasi-2D quantum well separated from 
the metal by a blocking barrier of thickness d. We as- 
sume that the quantum well and the barrier are char- 
acterized by the same dielectric constant e. We focus 
on the properties that are essentially independent of the 
specific atomic and crystal structure of the materials. 
Therefore, we assume the idealized model of the metal, 
which consists of a positive background charge and an 
arbitrary number of nearly free electrons, and its surface 
is the infinite plane. The single electron generates the 



Coulomb field, which induces the positive charge on the 
metal surface, which in turn acts on the electron chang- 
ing its quantum state. For our model structure the in- 
teraction of the electron with the induced charge can be 
calculated by the image charge method. We assume that 
the electron is strongly confined within the quantum well 
in the normal-to-plane (vertical) direction z and the im- 
age charge distribution is a mirror reflection of the elec- 
tron charge distribution with respect to the metal surface 
plane (z = 0). 

The potential energy of the interaction of the electron 
with the charge induced on the metal surface can be ex- 
pressed as 



W = 



d^rd 



,3 / e-(I■)^?+(r') 



(1) 



where prefactor 1/2 accounts for the self- interaction char- 
acter of the electron-image charge interaction and k — 
l/47r£o£- In Eq- O, there appear the two spatially 
separated charge densities: the electron charge density 
p_(r) — —e\tp(x,y)\'^S{z — d) and the image charge den- 
sity q+{y) = q\ip{x,y)\^6{z + d), where e > is the ele- 
mentary charge, _!/^a;, y) is the wave function of the elec- 
tron, and q = e 

First, in the framework of the simple 2D model, we 
perform the variational calculation with the trial wave 
function chosen in the form 



^{x, y) - {l/Ti^'^l) cxp[-(a;2 + y^)l2f 



(2) 



where variational parameter I determines the radius of 
electron localization. For the Gaussian charge densities 
we can apply the effective interaction potential ;19j when 
calculating the potential energy [Eq. d)]. Expression 
takes on the form (cf. Eq. (19) in Ref. '19'|) 

W{1) = -(l/2)(Ke(7//)(V2)'/' erfcx(2i/2d/Z) . (3) 

Using wave function (jJJ we calculate the expectation 
value of the total energy E{1) = h^/{2mj'^) + W{1), 
where rrie is the effective electron mass. Function E[l), 
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displayed in Fig. 1 for d = an, possesses the pro- 
nounced minimum for finite radius of electron local- 
ization / = 3.61au with E^iri = — 0.158i?£), where 
R]j — rrif.K^e^ / {2h?) is the effective donor rydberg and 
au = / (rrieKe^) is the effective Bohr radius. The re- 
sults of Fig. 1 show the existence of the bound state with 
the lateral electron localization. The electron moves par- 
allel to the interface not as a plane wave but as a wave 
packet of finite spatial extension in lateral directions. In 
order to conveniently picture this motion, we find it use- 
ful to introduce a new quasi-particle, which we call in- 
ducton, since the confinement potential, which forms this 
wave packet, is generated by the induced charge. 
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energies in the self-consistent method, e.g., Hartree-Fock 
method. We note that one can obtain the same Hartree 
potential Q using the quantum linear response theory. 
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FIG. 1: Expectation value E of the ground-state energy of the 
inducton calculated with trial wave function (2) as a function 
of radius I. Inset: Exact (dots) and variational Gaussian 
(solid curve, Eq. ^) ground-state wave function as a function 
of lateral position ry. 

For the stationary states the Hamiltonian of the in- 
ducton can be derived from the condition 5{{iI}\T\iIj) + 
W)/5'4>* — 0, where T is the kinetic energy operator and 
W = W['ijj* ,ip] is given by (^3). We obtain the Hamilto- 
nian of the form H = — ?i^V^/(2me) + U{r\^,d), where 

[/(r,| ,d) = -.eqj d^ (4) 

and r|| = (x,y). Let us discuss the physical content of 
Hamiltonian H. Quantity U [Eq. is the Hartree po- 
tential energy due to the charge induced on the metal 
(dielectric) surface, but is not the potential energy of the 
electron, which absolute value is two times smaller [cf. 
Eq. 1^]. As a result, eigenfunctions V'n of Hamiltonian 
H can be identified with the electron wave functions in 
the inducton states, but the eigenvalues of H, En, do 
not determine the total energy of the system. Instead 
En possess the same interpretation as the one-electron 
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FIG. 2; Ground-state energy: total (solid curve) and one- 
electron (dashed curve), and radius / (dotted curve) of the 
inducton as functions of distance d from the quantum well to 
metal surface. 

Due to the interaction with the induced charge the 
electron is not an isolated particle. During the transition 
between quantum states the induced charge can either 
transfer its energy to the electron (the transition with 
full relaxation) or release it partly or fully into the metal 
(the transitions with partial relaxation or without relax- 
ation). In the first case, we can observe the total energy 
of the inducton, while in the last case, we can measure 
the one-particle energy, similarly as in the many-electron 
system described by the Hartree-Fock method under con- 
dition that the Koopmans theorem is satisfied. Figure 2 
displays the one-electron lowest energy level, total energy, 
and radius I of the inducton. We note that the inducton 
radius is comparable with d. 

The results depicted in Figs. 1 and 2 have been ob- 
tained with Gaussian trial wave function In order 
to check the validity of our choice of wave function ^ , 
we have solved the eigenequation of H by the numer- 
ical self-consistent method jl^ on the one-dimensional 
radial mesh. In the present paper, we limit our atten- 
tion to the s-like solutions with the rotational symmetry. 
The numerical procedure yields the results, which can 
be treated as " exact" . We have found (cf. inset of Fig. 
1) that the ground-state wave function obtained by the 
self-consistent procedure is almost the same as the op- 
timum trial wave function [Eq. (2)], which means that 
the exact ground-state wave function can be very well 
approximated by the Gaussian of form This proves 
the consistency of our simple model based on Eq. |(5J). 

The electron motion parallel to the interface can be 
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described by a plane wave, but then the induced charge 
density is zero, total energy E > 0, and the inducton 
states are not formed. In the inducton states with E < 0, 
in which the electron wave function is localized, potential 
U [Eq. 10}] does not vanish and is attractive. In these 
states, the momentum of the electron is not defined. The 
inducton can move in the lateral direction with the non- 
zero momentum but the electron wave function in the 
inducton state is a wave packet. 

The inductons considered so far were formed in the 
quantum well; so, they were described as the 2D systems. 
Let us generalize our approach to the 3D systems. For 
this purpose we abandon the assumption of the quantum- 
well confinement and consider the electron moving in the 
half-space outside the blocking barrier. The electron in- 
duces the charge on the metal surface and due to the 
interaction with the induced charge the 3D wave packet 
is formed. The system is described by the correspond- 
ing 3D Hamiltonian, which has the form of its 2D coun- 
terpart with potential energy U{r) = C/(r||,z) resulting 
from the charge induced on the metal surface and poten- 
tial of the blocking barrier, i.e., U{r) = oo for < z < d, 
and U{r) is given by the 3D version of formula Q for 
z > d. The image charge density is the mirror reficction 
of the electron charge density in the z = plane, i.e., 
Q+A^) = lHn{x,y,-z)\^. 
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FIG. 3: Wave functions of the three lowest-energy s-like states 
plotted as functions of cylindrical coordinates ry and z. 

Assuming the cylindrical symmetry, we solve the corre- 
sponding 3D eigenproblem on the 2D mesh (r|| , z) by the 
self-consistent procedure. Figure 3 depicts the s-like wave 
functions, which correspond to the lowest-energy levels. 
The ground-state wave function (n = 1), as a function 
of rii , has the Gaussian shape. The excited-state wave 
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functions (n 
of variables r 



= 2,3) clearly exhibit the non-separability 
and z. 

Figure 4 shows that with the increasing blocking bar- 
rier thickness the energy levels shift upwards and the in- 
terlevel separations decrease. This feature is in a qual- 
itative agreement with the experimental data 0. We 
note that the energy levels do not form the hydrogenlike 
Rydberg series. 



We expect that the quasi-2D inductons can occur in 
planar semiconductor nanostructures, which contain the 
quantum-well layer and are covered by the metal layer. 
The largest binding energy of the inducton should be 
observed in semiconductors with small dielectric constant 
and large electron band mass. Due to the larger valence 
band mass, the hole inducton should be more stable than 
the electron inducton. 

The 3D inductons can appear in semiconductor struc- 
tures of the type of the field effect transistor with the 
isolated gate. Due to the increase of the binding energy 
with the decreasing blocking barrier thickness, the induc- 
tons formed in the structures containing the thin blocking 
barrier will be the most strongly bound. A possible ex- 
perimental evidence of such inducton states is within the 
reach of modern experimental setups. 
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FIG. 4: One-electron s-like energy levels calculated with e = 1 
as functions of thickness d of the blocking layer. Inset: Fourier 
transform of the ground-state inducton wave function (solid 
curve) vs fcy (in A~^). Dots correspond to the experimental 
data taken from Fig. 4(A) of Ref. Q for time delay 266 fs. 



However, the so-called image-potential electron states 
[lo| were observed |^ H, |1] near the metal-vacuum in- 
terface. The theoretical interpretation of these states is 
based upon the notion of Rydberg states 0, ■ In this 
interpretation 0, j the potential energy of the interac- 
tion of the electron with its image charge is assumed in 
the form Uimage{z) = —Keq/{A:z), which is valid for the 
classical point charge. In quantum mechanics, the charge 
density of the electron is determined by a wave function. 
The image of the spread-out charge is also spread out. 
The potential energy of the electron should be calculated 
according to Eq. |^ and the induced potential should 
depend on lateral position r\\ [cf. Eq. Q)]. Potential 
U image, uscd in Rcfs. 0, ll| , corrcspouds to either the 
classical particle or the quantum particle with the lateral 
wave function given by a Dirac delta function. In spite of 
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this the authors [7], Ig use the plane wave for the lateral 
motion of the electron, which is a serious inconsistency. 
If the electron travels in the x — y plane as the plane wave 
(cf. Fig. 1 for I — > oo), its image charge is completely 
delocalized and the induced potential is zero 0|. This 
inconsistency is repeated in all the other papers on this 
subject, e.g., i,[l3,[l3. 

Based on the results of the present paper, we propose 
an alternative interpretation of the states observed in 
Refs. Bli. According to our proposition these find- 
ings [l,!^] can be treated as the experimental evidence 
for the formation of the inducton. To support this inter- 
pretation we have calculated the Fourier transform of the 
ground-state wave function of the inducton. The results 
(inset of Fig. 4) very well agree with the experimental 
data 1^. Using the results of Fig. 4 for = 0, we have 
calculated the energy differences A^n — Em — En be- 
tween the one-particle energy levels and compared them 
with the experimental data 0. The calculated (mea- 
sured 13) energy differences are A21 = 0.41 (0.43) eV, 
A31 = 0.49 (0.50) eV, and A41 = 0.54 (0.56) eV (en- 
ergy separation A21 estimated from Ref. j5| is 0.45 eV). 
We see that the agreement with experiment is very good. 
The fact that the one-particle energy levels are observed 
in experiments 0, 0| can result from the slow relaxation 
of the induced charge during the photoemission. We note 
that the present model is valid if the electron does not 
penetrate the metal j^^, i.e., if the metal surface is cov- 
ered by a thin dielectric layer. This condition is satisfied 
in the experiments 

In summary, we have shown that the charge induced on 
the metal (dielectric) surface by the nearby electron cre- 
ates the lateral confinement potential, which in turn leads 
to the lateral localization of the electron. The collective 
electron-induced-charge motion can be conveniently de- 
scribed with the help of the new quasi-particle, the in- 
ducton. According to the reinterpretation we have pro- 
posed the localized states observed near the metal surface 
are the inducton states. The electron described by the 
plane wave does not generate any induced charge poten- 
tial. The plane-wave states can only be created if the 
excitation energy is larger than the continuum-energy 
threshold. The electron localization in the vertical di- 
rection is only possible if the lateral electron localization 
occurs. Electrons states which are not laterally localized 



cannot be bound to the metal surface. 

We expect the inducton states to be universal for 
the charge carriers moving near the metal or dielectric 
surfaces and should observed at the metal-vacuum and 
dielectric- vacuum interfaces as well as in semiconductor 
nanostructures. 

This work has been supported in part by the Polish 
Government Scientific Research Committee (KBN). 
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